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1 Problem 

A variant on the electro- or magnetostatic boundary value problem arises in accelerator 
physics, where a specified field, say B(0,0,z), is desired along the z axis. In general there 
exist static fields B(x,?/, z) that reduce to the desired field on the axis, but the "boundary 
condition" B(0, 0, z) is not sufficient to insure a unique solution. 
For example, find a field B(x, y, z) that reduces to 

B(0, 0, z) = Bo cos /czx + Bq sin kzy (1) 

on the z axis. In this, the magnetic field rotates around the z axis as z advances. 

The use of rectangular or cylindrical coordinates leads "naturally" to different forms for 
B. One 3-dimensional field extension of (|]) is the so-called helical wiggler |]], which obeys 
the auxiliary requirement that the field at z + 5 be the same as the field at z, but rotated 
by angle k5. 

2 Solution 

2.1 Solution in Rectangular Coordinates 

We first seek a solution in rectangular coordinates, and expect that separation of variables 
will apply. Thus, we consider the form 

B x = f(x)g(y) cos kz, (2) 

B x = F(x)G(y)smkz, (3) 

B z = A(x)B(y)C(z). (4) 

Then 

V • B = = f'g cos kx + FG' sin kx + ABC' , (5) 

where the ' indicates differentiation of a function with respect to its argument. Equation (|5|) 
can be integrated to give 

f 9 FG' 
ABC = — — sin kz-\ ; — cos kx. (6) 

k k 

The z component of V x B = tells us that 

-77— — fg' cos kz = ^—^ = F'G sin kz, (7) 
oy ox 
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which implies that g and F are constant, say 1. Likewise, 

— — = —fk sin kz = ^—^ = A'BC = — ^— sin kz, (8) 

oz ox k 

using (§-0). Thus, /" - k 2 f = 0, so 

/ = he kx + f 2 e~ kx . (9) 

Finally, 

^ = GA;cosA;z = ^ = y45'C=^-sin^, (10) 

SO 

G = Gxe ky + G 2 e~ hy . (11) 
The "boundary conditions" /(0) = B = G(0) are satisfied by 

f = B cosh fca;, G = B cosh fcy, (12) 

which together with (||) leads to the solution 

S x = i? cosh kx cos fcz, (13) 
5 y = B cosh fcy sin kz, (14) 
.Bz = — Bq sinh /cxsin /cz + B sinh fcy cos fcz, (15) 

This satisfies the last "boundary condition" that B z (0,0,z) = 0. 
However, this solution does not have helical symmetry. 

2.2 Solution in Cylindrical Coordinates 

Suppose instead, we look for a solution in cylindrical coordinates (r, 6, z). We again expect 
separation of variables, but we seek to enforce the helical symmetry that the field at z + 5 
be the same as the field at z, but rotated by angle k8. This symmetry implies that the 
argument kz should be replaced by kz — 6, and that the field has no other 9 dependence. 
We begin constructing our solution with the hypothesis that 

B r = F(r)cos(kz-6), (16) 
B e = G(r)sm(kz-9). (17) 

To satisfy the condition ([!]) on the z axis, we first transform this to rectangular components, 

B z = F(r) cos(kz - 6) cos # + G(r) sin(£;z - 0) sin 0, (18) 
By = ~F(r)cos(kz-0)sin0 + G(r)sm(kz-0)cos0, (19) 

from which we learn that the "boundary conditions" on F and G are 

F(0) = G(0) = B . (20) 
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A suitable form for B z can be obtained from (V x B) r = 0: 

1 dB z dB ft 



r 89 dz 



kG cos(kz — 9), 



so 



B z = —krG sin(kz — 9), 

which vanishes on the z axis as desired. 

From either (V x B) e = or (V x ~B) Z = we find that 

p d{rG) 



dr 



Then, V • B = leads to 



F = C- 



dh 



G = C- 

kr 2 



(krf 
1 + + 



' 3(kr) 2 
1 + - A rr Z - + 



(21) 



'.22) 



<*> w +fcr w- [1+(ir)21(fcrGH0 ' 

This is the differential equation for the modified Bessel function of order 1 [[|. Hence, 

t h{kr) C 



(23) 



(24) 



d(kr) V u At 7 2 
The "boundary conditions" (|20| ) require that C = 2£? , so our second solution is 



(25) 
(26) 



B r = 2B flo(fcr) - cos(^ 

B = 2S ^sin(A;2 - 9), 
kr 

B 7 = 



-2B I lS m{kz-9), 



(27) 

(28) 
(29) 



which is the form discussed in 
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